Dilute helium solutions in particle-physics experiments: response to a heat source 
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Motivated by a proposed experimental search for the electric dipole moment of the neutron 
(nEDM) utilizing neutron- 3 He capture in a dilute solution of 3 He in superfluid 4 He, we reana- 
lyze a measurement in a related system by Lamoreaux et al. of the 3 He density distribution in a 
dilute solution in the presence of a heat flow. We show, from considering the thermodynamics and 
hydrodynamics of dilute solutions, that this experiment in fact measures the thermal conductivity 
of the solution. Calculating the thermal conductivity within a microscopic framework, taking into 
account phonon-phonon, phonon- 3 He, and 3 He- 3 He scatterings, we find satisfactory agreement of 
theory with the data, serving to confirm our understanding of the microscopies of the helium in the 
nEDM experiment. 

PACS numbers: 67.60.G- 13.40.Em 
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I. INTRODUCTION 

The proposed experiment [l[ to measure the neutron 
electric dipole moment at the Oak Ridge National Labo- 
ratory Spallation Neutron Source (SNS) utilizes the ab- 
sorption of ultracold polarized neutrons on polarized 3 He 
atoms in solution in superfluid 3 Hc, via the reaction 



He -> p + t + 764 kcV. 



(1) 



A crucial issue is to be able to periodically sweep out the 
3 He by means of heat pulses |2J ; the underlying physics 
of the transport is scattering of phonons in the superfluid 
He against the 3 He. A pioneering experiment by Lam- 
oreaux et al. Q (hereinafter referred to as I), measuring 
the effect of a heat source on the steady state distribu- 
tion of 3 He atoms in a dilute solution in superfluid 4 He, 
provides a prototype for the effects of a phonon wind on 
the 3 He. 

The experiment described in I utilizes a novel tech- 
nique to measure the thermal transport properties of di- 
lute solutions. The solution is contained in a cylindrical 
cell roughly 5 cm in diameter and 5 cm long, cooled at 
one end by a dilution refrigerator to a temperature T in 
the range 0.45 — 0.95 K. The concentration of 3 He in solu- 
tion, X3 = 713/(71,3 + 714), where 773 and 774 are the 3 He and 
He number densities, is in the range 7 x 10~ 5 to 3 x 10~ 4 . 
A temperature gradient is created by a resistive heater 
generating 7-15 mW, located roughly midway between 
the ends and near the cylinder wall (see Fig. @] below). 
The resulting spatial distribution of the 3 He density is 
probed by a well-collimatcd neutron beam of diameter ~ 
0.25 cm going through the cell. A fraction of the neutrons 
are captured via the reaction ([1} . The XUV scintillation 
light from protons and tritons in the liquid 4 He is con- 
verted to visible light and detected by a photomultiplicr 
tube which views the solution through a window at the 
other end of the cell. The yield of scintillation light, 



measured as a function of cell temperature, initial 3 He 
concentration, and heater power is used to determine the 
thermally induced change in the 3 He distribution. 

The experiment was originally analyzed in terms of 
diffusion of the 3 He at fixed temperature, described by 
a diffusion constant. In this paper we reconsider the re- 
sults of the experiment at low temperatures in terms of 
the full hydrodynamical description of dilute solutions of 
3 He and superfluid 4 He In the temperature regime, 
T < 0.6 K, on which we focus, phonons are the dominant 
excitations in the superfluid, and carry most of the heat. 
In the following section we lay out the basic thermody- 
namics and hydrodynamics of dilute solutions, and show 
that under stationary conditions in the relevant range of 
£3, gradients of the 3 He density are, in fact, accompanied 
by non-negligible temperature gradients. We conclude 
that the experiment measures the overall thermal con- 
ductivity of the solution, rather than the 3 Hc diffusion 
constant. 

The basic scattering mechanisms determining the 
transport properties of dilute solutions - 3 He- 3 He, 3 He- 
phonon, and phonon-phonon interactions - are spelled 
out in Sec. III. In the concentration regime of the experi- 
ment in I, the 3 Hc- 3 Hc interactions are sufficiently strong 
that they keep the 3 He in thermal equilibrium at rest at 
the local temperature T(r). The phonon-phonon interac- 
tions keep the phonon momentum distribution in equilib- 
rium, drifting at a velocity v p h- In Sec. IV, we calculate 
the thermal conductivity (dominated by the phonons) 
in this situation, by considering the drag force on the 
phonons due to their scattering against the stationary 
3 Hc. The method we use is identical to that used in an 
earlier calculation of the mobility of ions in superfluid 
4 He and later used by Bowley in his accounting 
of the Lamoreaux et al. experiment. Our calculation of 
the momentum transfer between phonons and 3 He agrees 
in structure with Bowley's calculation. In Sec. V, we re- 
late the effective diffusion constant extracted in I to the 



2 



actual thermal conductivity of the dilute solutions, and 
show how the microscopic theory satisfactorily explains 
the experimental findings in I. 

At lower 3 He concentrations, x < 10~ 6 , the transport 
must be calculated by solving the coupled Boltzmann 
equations for the phonons and 3 He; the results will be 
given in Q. Finally, the impact of these results on the 
transport of 3 He in the much lower concentration regime 
of the SNS neutron electric dipole moment experiment 
will be described separately in a third paper 0] . 



II. RESPONSE TO STATIC DISTURBANCES 

To understand how dilute solutions of 3 He in super- 
fluid 4 He respond to a localized static heat source, we 
first review the hydrodynamics of the solutions. At low 
temperatures the phonons are the dominant excitations 
of the 4 Hc, and the momentum density or mass current 
density of the 4 He is 



.94 = P S V S + PphVph, 



(2) 



where v s is the superfluid flow velocity, v p h is the phonon 
fluid flow velocity, p s is the superfluid mass density, and 
Pp?i( < C Ps) is the 4 He normal fluid density. Similarly, the 
3 He momentum density is 



.93 = P3«3, 



where v$ is the 3 He flow velocity, and P3 



(3) 

i*na, with 

to* = TO3 + 5m ~ 2.34 TO3 the effective mass [10| . The 
total mass current, g, is 174 + 33. 

When the 4 He mass flow vanishes, \v s \ = 
(p p h/ Ps) v p h "C v p h at temperatures < 0.6 K, and thus 
in the absence of 3 He mass flow, the only relevant flow 
velocity is that of the phonons. Then force balance in 
the dilute solutions implies that to linear order, 



VP = T) p hV'Vph 



(4) 



where P = P3 + P p h is the total pressure, with P3 = n%T 
the 3 Hc partial pressure (we generally work in units with 
h and Boltzmann's constant, ks, equal to unity), T is 
the temperature, P p h is the phonon partial pressure, and 
rjph the first viscosity of the normal fluid. In the exper- 
iment in I, the viscosity term is insignificant compared 
with the drag of the 3 He on the phonons, as we show 
in Sec. HVl Thus the the total pressure is effectively con- 
stant throughout the system. (However, at the ultralow 
3 He concentrations of the proposed SNS experiment, the 
phonon viscosity does play a significant role |{|.) 

In addition, as one sees from the linearized superfluid 
acceleration equation Q , 



at 



0. 



(5) 



the 4 Hc chemical potential, fj,±, is also constant in a 
steady state. The Gibbs-Duhcm relation for the solu- 
tions, VP = n 4 V^ 4 + n 3 V^ 3 + SWT, with ^3 the 3 He 



chemical potential, T the temperature, and S the total 
entropy density, then implies that 



VP 3 + SphVT = 0, 



(6) 



where S p h = 4P p / l /P is the phonon entropy density, and 
dPph = SphdT. Equation ^ then tells us the simple 
relation between the temperature and 3 He density gradi- 
ents 



VP 



T 



Sph + m 



Vn 3 ; 



(7) 



in a steady state a gradient of the 3 He density is always 
accompanied by a gradient of the temperature. 

The heat flux, Q, in I, is related to the temperature 
gradient by Q = —KVT, where K is the thermal con- 
ductivity of the solution. Thus Eq. ([7]) relates the 3 He 
density gradient to the heat flux by 



„ S ph + n 3 -. 
3 = ~~TK — 



(8) 



We see that the basic transport coefficient measured in I 
is the thermal conductivity. A second feature is the role 
of the 3 He density. Since 



~>ph 



2tt^ (T 
45 



one has 



Sph P 



(9) 



(10) 



with T measured in K; at T = 0.45 K and i 3 = 3x 10~ 4 
the ratio is unity. 



III. MICROSCOPIC SCATTERING PROCESSES 



The microscopic processes that determine the trans- 
port properties of dilute solutions are 3 He- 3 He llj . 
phonon-phonon [T^], and phonon- 3 He [T^j scatterings; 



the amplitudes of these processes are all well known from 
earlier helium research. At low energies, the total cross 
section for a 3 He atom scattering from a second 3 He atom 
of opposite spin can be written as 



9tt 3 

CT 33 = ~T2~ V 
D 



33.0 



/»4 



/ TO4S 



D 



10.3A 2 , (II) 



where ko is the Debye wave number, defined by ?i 4 = 
k^l&TT 2 , and ^33,0 = —0.064 measures the strength of 
the 3 He- 3 He interaction at zero momentum transfer fllj ) . 
Phonon-phonon scattering rates are rapid compared with 
phonon- 3 He rates. While such scatterings do not directly 
affect the heat current, they do play the important role of 



3 



keeping the phonon momentum distribution, n#, in local 
thermodynamic equilibrium with a drift velocity, {f p /j, 



1 



,(sq-q-v ph )/T _ I ' 



(12) 



The amplitude for scattering of a phonon of wavevector 
q against a 3 Hc of wavevector p to final states q' and 
p' is determined theoretically in terms of the measured 
excess volume, a, of a 3 He atom compared with that of 
a 4 He atom and the 3 He effective mass, m* . To a first 
approximation the scattering is elastic, in the sense that 
q = q' and p = p', with amplitude 



(p'q'\T\ P q) = ^-(A + Bcos 



(13) 



where 9 is the angle between q' and q, A = n^da/dn^ 
and B = (l + a + 5m/m4 : )(m4 : /m*)(l + a — m 3 /mA). The 
momentum dependent scattering rate of phonons from 
the 3 He is then 



where 



73 (<?) = "j J, 



J = A 2 + B 2 /3~2AB/3; 



(14) 



(15) 



the q in this rate is characteristic of Rayleigh scattering. 

We estimate the importance of 3-3 versus phonon scat- 
tering in bringing the 3 He into equilibrium. The mean 
free path of a 3 He scattering on unpolarizcd 3 Hc is 



*33 = 



n 3 a 33 



8.66 x 10- 

%3 



(16) 



the factor n 3 / '2 is the density of opposite spin 3 He. Simi- 
larly, the mean free path for scattering of 3 He of momen- 
tum p on phonons is given by p 3 /mT, where 



r = -1 



d 3 q 

(2ny 



-q 2 l3 (q) n° (1 + n°) 



270ttJ| ^r-T 9 , 

7l4 J S Z 



(17) 



l) 1 is the 



in the limit p q. Here n° q = [e sq l T 
equilibrium phonon distribution function. Replacing p 2 
by 3m*T, an appropriate thermal average, we find 



t 3 ph 



Vs fs, 



2. J 



0.077 



>ph 

ri4 



T 3 / 2 



0.45 K 

T 



15/2 



(18) 



Comparing the mean free paths (|16p and (fl"8"|). we find 



l 3 ph 

£33 



0.89 x I0 b x 3 



0.45 if 
T 



15/2 



(19) 



For T = 0.45K and x 3 = 10" 6 , £ 3ph « £33, while for 
T = 0.65K and x 3 = 3 x 10~ 4 , t 3ph /i 33 « 16.9. Under 
the conditions of the experiment of I, a good first approx- 
imation is to assume that scattering of 3 He by 3 He atoms 
is fast compared with scattering of 3 He by phonons. In 
fact, as will be shown in detail in we may take the 
momentum distribution of the 3 Hc to be simply that of 
a classical gas in equilibrium at rest, 



f° 



-(p 2 /2m*-M3)/T 



(20) 



IV. THERMAL CONDUCTIVITY 



The thermal conductivity of the dilute solutions can 
in general be calculated by solving the coupled phonon 
and 3 He Boltzmann equations to determine the steady 
state momentum distributions of the phonons and the 
3 He, and, from these distributions, the heat and particle 
currents. Since the phonon heat current density, Q p h, is 
nothing but the phonon momentum density times s 2 , 

Q P h — s 2 [ ,t q ^ q\ 



(2n} 



^phVph t 



(21) 



where nq is the distribution of the phonons, the phonon 
thermal conductivity can be found simply under the con- 
ditions of the experiment in I by calculating the rate at 
which the phonons lose momentum by scattering against 
the stationary 3 He (f], Q . 

In a steady state the net force density on the phonons 
drifting at velocity v p h is balanced by the phonon pres- 
sure gradient, VP p h = S p hVT. Microscopically then, 



VPph 



d 3 q f d 3 q' 



d 3 p 



(2tt) 3 J (2tt) 3 J (2tt) 3 
xq\(T)\ 2 2TrS(sq + p 2 /2m* - sq' -p' 2 /2m*) 
x [/°n f (1 + n r ) - f p ,n r (1 + n*)] , 



(22) 



where p' — p = q — q' = k, (T) = (p'q'\T\pq), and the 
factor of 2 is from the 3 He spin sum. In terms of the 3 He 
structure function, 



S 3 (k,uj) 



d 3 p 



= n 3 



5(u+p 2 /2m* 
1/2 



P n l2m*) /o, 



27rfc 2 T 



which obeys S 3 (k,~uj) 
Eq. (EH) as, 



,-m* (ui-k 2 /2m*) 2 /2k 2 T 

(23) 

J / T S 3 (k,u)), we can rewrite 



VP ph = - 



d 3 q f d 3 



(2tt) 3 J (2tt) 3 
xq \(T)\ 2 2ttS(lj -sq + sq')S 3 (k,u) 

ng-(l + n$') - nq->(l + n^)e _ " /T 



(24) 
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FIG. 1: (Color online) Thermal conductivity of a dilute solu- 
tion of 3 He in superfluid 4 He at T = 0.45 K. The approximate 
phonon thermal conductivity, Eq. (|27[) (dot-dashed line), in 
the limit in which the 3 He are at rest in equilibrium is shown 
together with the full calculation from Ref. [?J (dashed line); 
these results differ only by a few percent at the lowest concen- 
trations. The solid line shows the phonon thermal conductiv- 
ity calculated with phonon recoil (see text). 
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FIG. 2: Ratio of the contribution of the viscous term (see 
Eq. ffl) to the drag on the phonons due to scattering against 
the 3 He at a temperature of 0.45 K. The magnitude of the 
drag is given by the 3 He pressure gradient (see Eq. (|26|H as 
a function of 3 He concentration X3 (the concentrations mea- 
sured in I are X3 = 7 x 10 -5 and 3 x 10~ 4 ). The viscous term 
is calculated for laminar flow in a circular pipe of radius R — 
2.5 cm, using the 4 He viscosity measured by Greywall [lr3 |. 



Since n f /(l + n f ) = n ? (l + n^)e^-^ h ^ T , the final 
line in Eq. (|2~4"|) to first order in v p h is n°(l + n° q ,)(k ■ 
Vph)/T. Finally, symmetrizing the integrand in q and q' , 
and carrying out the angular averages, we have, 



/CO 
du2Ti8(uj - sq + sq')S 3 (k 7 uj). (25) 
-00 

Neglecting 3 He recoil is equivalent to setting S 3 (k, w) = 
n 3 S(uj). In this case, 

x2nS(sq — sq'), 

n 3 T . 



3T 



Vph 



8tt 5 x 3 J /T n 8 



45 714 V s 



' P h. (26) 



Since VP p h = S ph ^T, and Q ph = TS ph v ph 
we conclude that 



Kph = 



ph 



n 3 T 90ttx 3 JT 



-K ph VT, 



(27) 



Figure [T] shows the phonon thermal conductivity, using 
the parameters A = -1.2 ± 0.2 [3, ^ , and B = 0.70 ± 
0.035 [l3|-[l5j]. The largest uncertainty is in A, owin g to 
a systematic difference between the measurements [13L 
[l4T | of the pressure dependence of the density of dilute 
solutions. Numerically, J = 2.2 ± 0.6. 

As Bowley emphasized @ , allowing for the recoil of the 
3 Hc in scattering produces a significant correction to the 
effective phonon- 3 He scattering rate (l~7| . In Appendix 
A we describe the calculations in detail. These effects 



increase the thermal conductivity by ~ 25 — 35% in the 
range T = 0.45 - 0.65 K (see Fig [J). 

We return now to justify neglecting the phonon viscous 
stress in the pressure equation ([?]). The crucial question 
is the magnitude of this stress on the phonons compared 
with the drag they experience from scattering against the 
3 He, Eq. (f2"5|) . We estimate the viscous stress, neglecting 
the drag, for laminar flow in a circular pipe of radius i?, 
where V 2 v p h is a constant, —8(v p h)/R 2 , throughout. As 
seen in Fig. the ratio is well below 1%. In fact, the 
viscosity for flow in a circular pipe plays a role only in 
a thin boundary layer, which for the concentrations and 
temperatures of the experiment in I is of order mm thick. 



V. COMPARISON WITH EXPERIMENT 

In the Lamoreaux et al. experiment, the heat flow pro- 
duces a non-uniform temperature and 3 He distribution, 
with the gradients of both T and n 3 proportional to the 
inverse of the total thermal conductivity of the dilute 
solution, Eqs. ([7]) and {5}. Note that the response of 
the 3 He is proportional as well to S p h + TI3. Since the 
3 He thermal conductivity is negligibly small throughout 
the range of concentrations in the experiment [owing to 
^33 *C i 3p h, Eq. (HHJ)], the experiment in fact measures 
the phonon thermal conductivity as limited by scatter- 
ing against the 3 He. In the calculation of the phonon 
thermal conductivity, deviations from the limit where 
3 He- 3 He scattering is dominant occur only at the low- 
est concentrations in the experiment and there only at a 
percent level. 

The response of the 3 He can be described in terms of 
the 3 He diffusion constant, D 3 , by noting that VP p h = 
— VP3, so that the relation between the phonon drift 



FIG. 3: (Color online) Contribution of the thermoelectric 
transport term relative to the effective diffusion constant, Dl 
extracted in I. The temperatures are 0.45 K (dotted), 0.55 K 
(dashed) and 0.65 K (solid). 



and the phonon pressure gradient can be written as (cf. 
Eq. (USD) 



Ih 
T 



VP 3 = n 3 v p h, 



(28) 



where in the limit that the 3 He are at rest in equilibrium, 



D, 



3T 2 



(29) 



In the experimental situation, the phonon drift drives 
VP3 rather than simply Wn 3 . Equation (|28p is, more 
generally, the response of the 3 He current to the phonon 
wind 0, 



:)3 = n 3 v 3 = n 3 v ph - D 3 Vn 3 - D T X7T, 



(30) 



in the situation where v 3 = 0. Here Dt = n 3 D 3 /T is an 
effective "thermoelectric" transport coefficient. 

The analysis in I, in contrast, neglects the temperature 
gradient, and determines an effective diffusion constant, 
Dl, by setting a simplified expression for the 3 He current 
to zero, 



j 3 L = n 3 v ph - D L \7n 3 = 0, 



(31) 



(cf. Eq. ([30]) ) thus defining Dl- Taking Q = TS p hV p h, 
Ref. I derives 



Vn 3 



n 3 



TS p hDl 



(32) 



To calculate Dl within our framework we compare 
Eq. 432) with Eq. © and find 



Di 



n 3 



K 



S p h + n 3 S p h 



(33) 



We can also write Dl in terms of the general 3 He trans- 
port coefficients as 



D, 



D 3 



Sph + n 3 ' 



-Di 



Sph 

Sph + n 3 ' 



-Da 



(34) 



b) 




FIG. 4: (Color online) Example of the 3 He distribution for 
a cell 5.2 cm in diameter and 5.0 cm high, with refrigera- 
tor temperature T r = 0.45 K, average £3 = 3 x 10 -4 , and 
heater power 15 mW. The sections are: a) perpendicular to 
the cylinder axis and containing the neutron beam axis; and 
b) perpendicular to the neutron beam (the cross in the figure) 
and through the center of the heat source. We assume in this 
simulation that the cylindrical wall of the cell is at the same 
temperature as that of the base where the refrigerator is at- 
tached; the top surface is insulated. The 3 He concentration 
contours, of constant spacing, are marked in units of 10~ 4 ; 
the minimum is £3 = 1.39 x 10~ 4 at the surface of the heater 
and the maximum is 3.03 x 10~ 4 on the lower cell boundary. 
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The Dt term in this equation makes a significant con- 
tribution at higher concentrations (see Fig. [3]); as noted, 
"3 ~ S p h for a significant part of the range of concentra- 
tions in I. 




■£ 0.85 
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FIG. 5: (Color online) Representative results of the finite 
element calculation of the relative integrated 3 He (column) 
density along the neutron beam as a function of the quantity 
£l = T 7 V I SphT where V is the heater power, along with data 
from Ref. 0]. The curves are: solid, 13 = 3 x 10~ 4 and fixed 
refrigerator temperature T r — 0.45 K with the temperature 
of the cell barrel, the "conducting" boundary, also fixed at 
0.45 K; dashed, the same calculation with only the refrigera- 
tor end at 0.45 K and the rest of the cell surface insulated; 
dotted, 13 = 7x 10~ 5 , T r — 0.45 K and the conducting bound- 
ary condition; the results for 13 = 3x 10~ 4 , T r = 0.65 K and 
the conducting boundary condition are, coincidentally, essen- 
tially the same as the dotted curve. The curves are plotted for 
the average temperature (greater than the boundary temper- 
ature) along the neutron beam path, resulting, at £l = 2.5, 
in positive shifts of £l ranging from A£l ~ 0.05 k 7 B cm 3 /s for 
the solid curve to A£ L « 0.5 k 7 B cm 3 /s for the T = 0.65 K 
case. 



An example of the 3 He distribution resulting from the 
finite element calculation is shown in Fig. UJ representa- 
tive results of integrated 3 He densities along the neutron 
beam are shown in Fig. [3] together with the data in Fig. 
4 of I. The solution depends quite sensitively on whether 
one considers the barrel of the cylinder to be thermally 
conducting or insulating; calculations with these two as- 
sumptions bracket the data. To illustrate the effect of 
other variables in the problem, we also show in Fig. [5] 
the results for T = 0.65 K and for the two concentra- 
tions used in I. 



VI. SUMMARY 

We have analyzed the heat current and the resulting 
3 He spatial distribution in dilute solutions of 3 He in su- 
pcrfluid 4 He in the temperature range where phonons are 
the significant excitations of the 4 He. We find that in the 
range of 3 He concentrations in the Lamoreaux et al. ex- 
periment 3, 7 x 10~ 5 < x 3 < 3 x 10~ 4 , heat transport 
is dominated by phonons. The response of the system is 
generally not a simple diffusion problem, described only 
by a 3 He diffusion constant, D3, since the temperature 
gradients are important. The experiment can be charac- 
terized as measuring the phonon thermal conductivity. A 
calculation of the 3 He density along the neutron beam in 
the geometry of the experiment indicates that the well- 
tested microscopic theory of 3 He-phonon scattering in di- 
lute solutions of 3 Hc in supcrfluid 4 He is consistent with 
these results as well. 



To further illustrate the physics, we now determine 
the 3 He density distribution for several limiting cases in 
the geometry of the experiment in I [l8j . There, the 
n.3 distribution was determined by solving V 2 u p h = 
assuming constant temperature throughout the volume. 
Here we determine the 77,3 distribution by first writing 



V • Q = -V • (K ph "VT) = -CV 



/VT 



= 0, (35) 



where from Eq. {57}, C 
P3 using the solution 



ls 2 /90n,J. We then eliminate 



SphT 



P, 



(36) 



of Eq. (|6]), where the constant P is the total pressure of 
the excitations. Equation (|35|) then reduces to a partial 
differential equation in T alone, which we solve using the 
finite element code FlexPDE. Given T(r) we then deter- 
mine n^(f) from Eq. (f3"6")> . and determine the constant P 
by fixing the total number of 3 He atoms in the system. 
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Appendix A: Recoil corrections 

In this Appendix we calculate contributions to the 
thermal conductivity due to the finite 3 He mass. We 
start from Eq. (|2"5"|) with the structure factor (|2"3"]) of the 
3 He, 
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VP, 



ph. 



_ 2Tm 3 f d 3 q f d 3 q' „ „ - 2 / m* \ 

~ v ph—QjT j y ^,n,(l + v)* K r )l ^^FtJ 

27rn 3 /■ d 3 q f d 3 q' k 2 \(T)\ 2 



6T J (2tt) 3 J (2tt) 3 4 sinh(sg/ 2T ) sinh(sg'/2T) 



1/2 



-m*(s(3r-sq'-fe 2 /2m*) 2 /2fe 2 T 



/ m* \ m * (;jg _ S9 ') 2 /2fc 2 T 
\2nk 2 TJ 



-fe 2 /8m*T 



(Al) 



r 



The integral in Eq. (|A1|) is proportional to the inverse 
of the thermal conductivity. In the limit m* — > oo 
the expression in square brackets reduces to S(sq — sq'). 
For finite m* there are two effects. First, the struc- 
ture factor, and therefore also the scattering rate, is re- 
duced in magnitude because of the nonzero momentum 
transfer, fe = q — if, as is shown by the final Gaussian 
factor. Second, as the first Gaussian factor indicates, 
there is an energy transfer sq — sq' which is of order 
(Tk 2 /™,*) 1 / 2 ~ (T/m*s 2 )^ 2 T. 

The contributions to the scattering amplitude for 
nonzero energy transfer and for nonzero velocity of the 



3 He atoms have not been investigated in detail, although 
the basic processes were discussed in Ref. [l(| ■ Here we 
make the simplest generalization of Eq. (|T3"|) that is sym- 
metric under exchange of q and g', 



s(aa') 1/2 

(p'q'\T\pq)= { ™ ) (A + Boost 
2n 4 



(A2) 



With prefactors omitted, the quantity to be calculated is 
thus 



J 



dq / dq' 



k(qq') 3 (A + Bcost 



1 hlnn 1 ^ 

dcosO — 

_ x 4smh(sq/2T)smh(sq'/2T) 



e -m' s 2 (q-q') 2 /2k 2 T e -k 2 /8m'T 



(A3) 



r 



We have evaluated the integrals numerically and find, as 
stated in Sec. IIV1 that inclusion of recoil effects increases 
the thermal conductivity by ~ 25 — 35% in the temper- 
ature range 0.45 — 0.65 K. 

The leading corrections to the result for m* — > oo are 
of order T/m*s 2 relative to the result in the low temper- 
ature limit; we now calculate them analytically. To first 
order in T/m*s 2 the effects of the nonzero momentum 
transfer and the nonzero energy transfer are additive, 
and we calculate each in turn. The more important term 
is due to the momentum transfer. When this term alone 
is included one finds 



Km(TK)/TK~ I 



1 (fc 4 ) 

8m*T (fc 2 ) 



where 



dq 



dcosO- 



BcosO) 2 



Asmh 2 {sq/2T) 



(A4) 



(...). (A5) 



In these integrals we may replace k 2 by its value 2q 2 (l — 
cos9) for zero energy transfer. The integrals over q and 
9 decouple and one finds 



TK 



T J I w 

4m* s 2 J I s 



lim (TK), 



(A6) 



where 



dcosO 



(A + Bcos6>) 2 (l-cos6>) 2 



U 2 Ub 



15 



:B 2 



and 



In = 



dq — - 



o 4sinlr (x/2) 



i\C(n), 



(A7) 



(A8) 



where £(n) is the Riemann zcta function of order n. 
Therefore the thermal conductivity is given by 



K 



25tt 2 J T 



1 



II J m*s 2 / Tm 



lim {TK). 



(A9) 



We now calculate the leading correction to the thermal 
conductivity due to the energy transfer, which is found 
by neglecting the term k 2 /8m*T in the exponent in Eq. 
(|A3|) . The Gaussian in the energy difference has a width 
small compared with T, so we adopt a procedure similar 
to that used in making the Sommerfeld expansion for 
low temperature Fermi systems, where the derivative of 
the Fermi function approaches a delta function. In an 
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integral of the form, 

roc e ~(x~y) 2 /2A 2 

G(x)=j o dy 9 (y) (2?rA)1/2 , (A10) 

where A is a constant, and the function g(y) varies slowly 
on the scale A, one finds on expanding g in a Taylor series 



about y = x, that 

G(x)=g(x) + ^-g"(x) + ... . (All) 

When this result is applied to the q 1 integral in Eq. (|A3[) 
with the final Gaussian omitted, one finds 



dq 



2q s 



dq / dq 
Jo 

T 



, f 1 dcosd ( m*s 2 \ 1/2 k 2 (qq') 3 (A + Bcos9) 2 e' m ' s2 ^') 2 /2k 2 T 



2 V 27rTfc 2 



4sinh 2 (sg/2T) 



4sinli (s<7/ 2r ) 2m* s 2 J 



dq 



1 dcos9 k 2 q 3 k 2 {A + B cos 9) 2 ( d 2 
-i 2 



k 2 q 



2smh(sq/2T) \d{q') 2 2smh(sq' /2T) J ql=q 



, (A12) 



which, when expressed in terms of the variables x = sq/T,y = sq' /T, and z = cos#, is proportional to 



J 



dx- 



2x* 



T 



dx — / — — 

'4sinh 2 (.T/2) m*s 2 J 2 sinh(a:/2) } _ x 2 

where n 2 = x 2 + y 2 — 2xyz. The second derivative is 

f 



{l-z){A + Bzf 



if 



2„,3 



K y 



dy 2 2sinHy/2)J y=x ' 



sinh(a;/2) 



1 + z) 



12 - jx 2 ~ Ax coth(.x/2) + -x 2 coth(x/2) 2 



(A13) 



(A14) 



On evaluating the integrals, one finds that the effect of 
energy transfer produces contributions to the integral 
[51 having the form 



(25tt 2 - 198) J 
33 



T 



(A15) 



The numerical factor (257r 2 — 198)/33 is approximately 
1.477 and thus one sees that the effects of nonzero en- 
ergy transfer are much less important than those due to 
nonzero momentum transfer. When both contributions 
to the thermal conductivity are included, one finds on 
inserting Eq. (f2"T|) for the low temperature limiting be- 
havior, 



K 



714s 



T 



90itx 3 JT 



IOOtt 2 - 198 J i 
33 J ~ 



(A16) 



The coefficient of J/J is approximately 23.90. For A = 
— 1.2 and B — 0.70, J is 2.16, J is 3.30 and their ratio is 
1.52. Thus the coefficient of T/m*s 2 is 35.5. Even though 
m* s 2 = 48.1 K, the effects of recoil are large even at 
temperatures well below 1 K as a consequence of the large 
numerical coefficient. This coefficient reflects the fact 
that the most important contributions to the momentum 
transfer arise from pho irons with momenta 3> T/s. At 
T = 0.5 K the correction is 37%, which is considerable. 
Higher-order contributions in T can be significant, and 
one would expect these to reduce the deviation from the 
low-temperature limiting result by an amount of relative 
order (0.37) 2 - 10%. These analytic results are m good 
agreement with the numerical integration of Eq. (|A3[) 
described above, and shown in Fig. [1] 
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